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Abstract: 

The accurate determination of the lowest electron attachment {EA) and ionization {IP) energies for molecules embedded in molec¬ 
ular junctions is important for correctly estimating, e.g., the magnitude of the currents (/) or the biases (V) where an I — V-curve 
exhibits a signihcant non-Ohmic behavior. Benchmark calculations for the lowest electron attachment and ionization energies 
of several typical molecules utilized to fabricate single-molecule junctions characterized by n-type conduction (4,4’-bipyridine, 
1,4-dicyanobenzene, and 4,4’-dicyano-1,1’-biphenyl) and p-type conduction (benzenedithiol, biphenyldithiol, hexanemonothiol, 
and hexanedithiol] based on the EOM-CCSD (equation-of-motion coupled-cluster singles and doubles) state-of-the-art method 
of quantum chemistry are presented. They indicate signihcant differences from the results obtained within current approaches 
to molecular transport. The present study emphasizes that, in addition to a reliable quantum chemical method, basis sets much 
better than the ubiquitous double-zeta set employed for transport calculations are needed. The latter is a particularly critical issue 
for correctly determining EA’s, which is impossible without including sufficient diffuse basis functions. The spatial distribution 
of the dominant molecular orbitals (MO’s) is another important issue, on which the present study draws attention, because it 
sensitively affects the MO-energy shifts <I> due to image charges formed in electrodes. The present results cannot substantiate the 
common assumption of a point-like MO midway between electrodes, which substantially affects the actual <t>-values. 
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1 Introduction 

Electron or hole injection into molecules embedded between 
two electrodes represents an important issue in the fabrica¬ 
tion of molecular devices. The efficiency of the charge in¬ 
jection and transport in molecular junctions is controlled by 
the highest occupied or lowest unoccupied molecular orbital 
(HOMO or LUMO, respectively), whichever is closest to the 
electrodes’ Eermi energy Ef, and the key quantity is the en¬ 
ergy offset eo = min(£'F - Ehomo,Elumo - Ef)- It can be 
compared to a tunneling barrier, which charge carriers have 
to overcome to generate a current. Ultraviolet photoelectron 
spectroscopy (UPS)—, thermopower— ■^, and transition voltage 
spectroscopy (TVS)—“—represent current methods to estimate 
the relative alignment of the dominant molecular orbital from 
experimental data. 

Recent analysis of a variety of transport data demonstrated 
that full current-voltage I — V curves beyond the ohmic regime 
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can be quantitatively reproduced by assuming that molecular 
transport is dominated by a single level (to be identihed with 
HOMO or LUMO) ^'^'*^~ — . This is an enormous simplihca- 
tion. Still, the correct description of the relative alignment Eq 
(= Em MO — Ef ot Ef — Efiomo) remains an important chal¬ 
lenge for ab initio approaches to the charge transport through 
single-molecule junctions. It is a challenge particularly be¬ 
cause of the high accuracy needed. Values directly determined 
in ultraviolet photoelectron spectroscopy (UPS) experiments 
amount to Co ^ 1 eV-. Results based on TVS by using a model 
able to excellently reproduce current-voltage (/ — V) curves 
measured in single-molecule junctions demonstrate that Co can 
be even smaller (Eq ~ 0.6 eV— ). It should be clear that, in 
view of such Iow-Eq values, estimates for /P’s and EA’s with 
errors ~0.5 eV typical for quantum chemical methods of mod¬ 
erate accuracy are unacceptable. Noteworthy, the quantity Eq 
is important not only because it determines the magnitude of 
the currents, but also because it indicates the biases beyond 
which an / — y curve becomes signihcantly nonlinear—. 

An important message, which the present study aims to 
convey, is that the accurate determination of the HOMO and 
LUMO energies (or, more precisely, the lowest ionization IP 
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and electroaffinity EA with reversed sign) represent a nontriv¬ 
ial issue even for isolated species of interest for molecular 
transport. Obviously, this is a minimal requirement for any 
molecular transport approach. 

In the present paper we report results of ongoing work fo¬ 
cusing on several prototypical molecular species, which are 
mostly utilized in the fabrication of single-molecule junctions, 
and examine the reliability of the results for ionization ener¬ 
gies IP{7ii —EffOMo) electroaffinities EA{f^ —Elumo) 
obtained within current methods employed in molecular trans¬ 
port with accurate estimates obtained within well-established 
quantum chemical methods. To avoid misunderstandings, let 
us explicitly mention that we will restrict ourselves through¬ 
out to the lowest IP and EA-, in fact, if at all, for the molecules 
analyzed below it is only the lowest EA that corresponds to a 
stable anion {EA > 0). 

The molecular species to be considered will include, be¬ 
sides 4,4’ -bipyridine (44BPY) — a molecule recently consid¬ 
ered from a TVS perspective —, molecules often used 

to fabricate two classes of molecular junctions. One class is 
represented by molecules characterized by n-type (LUMO- 
mediated) conduction, the other class comprises molecules 
characterized by p-type (HOMO-mediated) conduction. As 
specific examples belonging to the first class, we will consider 
44BPY1, BDCN (1 ,4-dicyanobenzene)i4 and 2BDCN (4,4’- 
dicyano-l,r-biphenyl). From the second class, we will ex¬ 
amine oligophenylene dithiols^ and alkanemono- and dithi¬ 
ol sii. 

Because of well-documented shortcomings of ubiquitous 
methods based on density functional theory (DFT), empha¬ 
sis will be on post-DFT methods. Besides methods al¬ 
ready utilized in approaches to molecular transport (e.g., 
A-SCF— , GW— and MP222), we will consider truly ab 
initio methods used in quantum chemistry: outer valence 
Green’s functions (OVGF) ^^'^^ , second-order algrabraic- 
diagrammatic constructions [ADC(2)] ^^4^ , and equation-of- 
motion (EOM) coupled-cluster (CC)^ approaches. The im¬ 
portance of using appropriate quantum chemical methods will 
be emphasized. As is well known, calculating EA’s is a very 
delicate problem in quantum chemistry altogether—. 

2 Methods 

Active molecules embedded in molecular junctions can be 
treated at various levels of theory ranging from tight-binding 
(extended Hiickel) to and post-DFT. In order to facilitate un¬ 
derstanding the message, which the results reported below aim 
to convey, we will briefly present the methods utilized in this 
study. 

(i) If the picture based on the self-consistent field (SCF) 
were valid (or, equivalently, electron correlations were ab¬ 
sent), the energy of highest occupied Hartree-Fock (HF) or¬ 


bital (i.e., HOMO) with reversed sign would represent the low¬ 
est ionization energy IP = —Escf,homo (Koopmans theorem). 
The physical meaning of the virtual (unoccupied) is contro- 
versia l^^i^^ . Virtual HF orbitals might have physical meaning 
if descriptions based on small basis sets succeeded (at least 
semi-)quantitatively, but this is often not the case. Table[^il- 
lustrates this failure, where results relevant in connection with 
existing MP2-based transport approaches and minimal (STO- 
6G) basis sets^ are presented. At the other extreme, it is also 
well known that the HF LUMO energy goes to zero in the 
complete basis set limi t^^i^* . This is a reason why attempts to 
“improve” the quality of a transport approach by using larger 
and larger basis sets end up with unphysical results^^. For 
large atomic orbital (AO) basis sets, the virtual HF orbitals 
have mathematical rather than physical meaning, namely, in 
providing an expansion manifold for the physical states of in¬ 
terest. In some cases, individual virtual or unoccupied HF 
orbitals (in particular, the LUMO) can reasonably describe, 
e.g., anionic bound or resonance states semi-quantitatively, 
provided that the size of the AO basis used is not too largei^. 

(ii) Within DFT-approaches, the single particle solutions 
of the Kohn-Sham (KS) equations are handled as if the 
corresponding eigenvalues/eigenfunctions were real orbital 
energies/wave functions. The corresponding implementation 
in a Landauer-NEGF formalism is straightforward because 
the DET description is mathematically a single-particle 
description. This is why the DET approaches to molecular 
transport are by far the most popular to date. Drawbacks 
of such DET-approaches are well documented. The drastic 
underestimation of the HOMO-LUMO gap and the related 
lineup problem (HOMO/LUMO energies too close to elec¬ 
trodes’ Eermi energy) are issues most frequently mentioned 
and not at all surprising: as is well known—, KS “orbitals” 
are mathematical objects rather than physical orbitals. 

The quantum chemical methods used in the present paper 
to compute the lowest electroafflnity (EA) and ionization IP 
energies are: 

(iii) The outer valence Green’s function (OVGF) 
metho d^^'^^ represents the most elaborate quantum chemical 
approach based on a single-particle picture. To the best of 
our knowledge, the OVEG approach has not yet been utilized 
in molecular transport studies. Therefore, let us mention 
that the OVGE method is a way to approximately include 
the contribution of the electron-electron interaction beyond 
HE. Details can be found in ref. . As explained below, 
this method is superior to the MP2-like approximation used 
recently— (see sec. 14.31 1. The OVGE method^2i22 exactly 
treats the full second- and third-order terms in the self-energy 
entering the Dyson equation for the one-electron Green 
function, and is augmented by a geometrical approximation 
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to also include further higher-order corrections^. Results 
for the electron affinities and ionization energies obtained by 
considering the second- and third-order terms are shown in 
the tables presented below (they are labeled as “2nd-order 
pole” and “3rd-order pole”, respectively) along with those of 
the full OVGF. The corresponding pole strengths^ are also 
indicated (percents in parentheses). 

(iv) The algebraic-diagrammatic construction (ADC) is 
based on a diagrammatic perturbation expansion. ADC(n) de¬ 
fines an approximate scheme of infinite partial summations 
exact up to the n-th order of perturbation theor y Its 

second-order version ADC(2) is superior to the so-called GW - 
approximation—, since only bubble contributions are included 
within the GW and not all second-order terms: expressing the 
self-energy E by the product of the single-electron Green’s 
function G and the effective interaction VT (E ^ GW, thence 
the name GW), vertex corrections^ are neglected within GW. 

(iv) Equation-of-motion (EOM) coupled-cluster (CC) ap¬ 
proaches at singles and doubles (CCSD)^, hybrid (CC2)— , 
and perturbative [CCSD(2)]^2i^ levels will be extensively 
applied in this paper. Corrections due to triples [EOM- 
CCSD(T)] will be also considered; the fact that these correc¬ 
tions are altogether negligible is an indication on the accuracy 
of the state-of-the-art EOM-CCSD method. 

(vi) Energy difference (A-) methods will be utilized for all 
the aforementioned cases. Within these methods, the low¬ 
est ionization energy IP and electron affinity EA are esti¬ 
mated as differences between the ground state energies S' of 
the corresponding molecular charge species (M=SCE, DET, 
MP2, CCSD, CC2) at the equilibrium geometry of the neu¬ 
tral molecule 

^ = '^M,cation ^M,neutral i (1) 

~ ^M,neutral ~ Sm ^anion ■ (2) 

3 Computational details 

The results of the SCE, DET/B3LYP, MP2, and OVGE cal¬ 
culations reported below were done with GAUSSIAN 09—. 
Coupled-cluster calculations of the IP, EA, total energies of 
the various charge species and excitation energies were per¬ 
formed with CEOUR— . Calculations within the so-called 
regular (strict) ADC(2) reported here have been done with 
the fully parallelized PRICD-E(2) code4^, which is inter¬ 
faced to MOLCAS— . As amply documented by extensive 
work of the Heidelberg theoretical chemistry group, the re¬ 
sults based on the strict ADC(2) are comparable to the second- 
order approximate coupled cluster singles and doubles model 
(CCiy^. Augmented with extra terms in an extended version 
[ADC(2)x]— , the results become comparable to the equation 
of motion coupled cluster singles and doubles method (EOM- 
CCSD)^. Unfortunately, a code enabling computations for 


molecular sizes of interest for molecular transport is not (yet) 
available. 

The inspection of the tables presented below reveals that 
the results obtained via the state-of-the-art IP- and EA-EOM- 
CCSD metho d and aug-cc-pVDZ (Dunning augmented 
correlation consistent double zeta) sets can be trusted. This 
is illustrated both by the good agreement between the EOM- 
CCSD and the A-CCSD values and by the fact that correc¬ 
tions due to triples [CCSD(T)] yield changes that are irrele¬ 
vant within numerical errors. Eor understanding the impact 
of polarization and diffuse functions we also present results 
obtained by using other basis sets: cc-pVDZ and cc-pVTZ, 
Pople basis sets (6-3IG*, 6-31IG* 6-31 l-H-G(d,p)),Dunning- 
Huzinaga double- and triple-zeta (DZ (DZ95 in ref.—), DZP, 
TZ2P), Karlsruhe basis sets (svp, dzp, tzp, qz2p), and (merely 
to compare with earlier MP2-based transport calculations, 
cf. Table© STO-6G. 

4 Results and discussion 

Because the tables presented in this paper, which contain very 
detailed information on both the methods and the basis sets 
employed, are self-explanatory, below we will only briefly 
emphasize the main aspects related to the lowest electron at¬ 
tachment (EA) and ionization (IP) energies of the molecules 
of interest. Still, as a technical remark, let us mention that, 
in view of the fact that double zeta (DZ) sets are ubiquitous 
in transport studies, among other basis sets, we have always 
included DZ-based results in the relevant tables. 

4.1 Lowest electron attachment energies 

Table [1] © and |3] collect results on the lowest electron at¬ 
tachment energies for 4,4’-bipyridine (CgH 4 N 2 , 44BPY), 1,4- 
dicyanobenzene (NC-CgH 4 -CN, BDCN), and 4,4’-dicyano- 
1,1’-biphenyl (NC-(CgH 4 ) 2 -CN, 2BDCN). 

Table[T]presents very detailed numerical results for 44BPY, 
a showcase molecule 13,15,46- ^^ jjj order to illustrate the main 
issues, which we have encountered in calculations of elec¬ 
troaffinities for molecules utilized in molecular electronics. 
The EA-EOM-CCSD method predicts a weakly bound anion 
44BPY*^ (EAZO). The essential condition for this is the in¬ 
clusion of a sufficient number of diffuse basis functions. As 
emphasized recently—, it is not the basis set size that mat¬ 
ters: as visible in Table [T] the basis set CC-pVTZ, which is 
larger than aug-cc-pVDZ, cannot stabilize the anion, just be¬ 
cause diffuse functions are missing. 

Notice also that a correct description requires a proper 
treatment of electron correlations. Even including sufficient 
diffuse basis functions (aug-cc-pVDZ), the SCE description 
(both Koopmans theorem and A-SCE) is even qualitatively in¬ 
adequate; the anion is predicted to be unstable. This is already 
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known from earlier worki^2il2iS. As visible in Table [T] we 
found that electron correlations (which, by definition, measure 
deviation from SCF) cannot be adequately included via MP2 
and the OVG F^^^^ : the 44BPY*^ anion remains unstable. 

On the other side, while agreeing among themselves, the 
CC2, ADC(2), and A-DFT all overestimate the anion stability. 
Concerning CC2, one can remark that A-CC2 performs better 
than EOM-CC2. 

Basically, the conclusions formulated above by analyzing 
the electron affinity of 44BPY also holds for the other two 
molecules [1,4-dicyanobenzene (BDCN) and 4,4’-dicyano- 

1.1 ’-biphenyl (2BDCN)] investigated; see Table |2] and [3 
Quantitatively, there is an important difference; the BDCN*^ 
and 2BDCN*^ anions are substantially more stable: their 
electroaffinities are by about 0.7 eV larger than that of the 
44BPY*^ anion. 

To summarize, the anions considered above can be accu¬ 
rately described at the CCSD level of theory provided that the 
basis set employed includes sufficient diffuse functions: triplet 
corrections [i.e., CCSD(T)] merely yield modifications of the 
EA-values within numerical errors, and the EA-EOM-CCSD 
values agree well with the A-CCSD values. 

4.2 Lowest ionization energies 

As representatives of molecules embedded in nanojunctions 
exhibiting a HOMO-mediated (p-type) conduction, we have 
studied and present detailed results for benzenedithiol (HS- 
CgI- 14 -SH, BDT, Table |4]i and related molecules (S-CgH 4 -S 
and S-C 6 H 3 E-S, Table [3), dibenzenedithiol (HS-(C 6 H 4 ) 2 -SH, 
2BDT, TablelS, 1,6-hexanemonothiol (H-(CH 2 ) 6 -SH, C6MT, 
TablelUl, and 1,6-hexanedithiol (HS-(CH 2 ) 6 -SH, C6DT, Table 
0 . 

As visible in these tables, the lowest ionization energies 
can be estimated with a good relative accuracy (^4%), which 
is satisfactory for quantum chemical calculations for many 
purposes, within IP-EOM-CCSD calculations by using rather 
modest basis sets. Although not dramatically large, the cor¬ 
responding absolute error (.$,0.4 eV) is still non-negligible 
from a molecular transport perspective in view of the rather 
small energy offset of the dominant molecular orbital relative 
to electrodes’ Eermi level. So, good basis sets are required 
not only for EA’s, but also for an adequate IP’s. Even for 
cations, the various lower level many-body approximations 
(MP2, CC2, ADC(2), CCSD(2) as well as their A-versions), 
deviating by up to ~ 0.4 eV from the EOM-CCSD approach, 
are still not too satisfactory. 

In view of the present results, the A-DET method cannot be 
recommended: deviations from the /P-EOM-CCSD estimate 
can be very large; the example presented in Table 0 indicates 
an error of ~ 0.8 eV. 

A special mention deserves the OVGE approximation. 


which appears to provide the /P-estimates closest to EOM- 
CCSD; the differences are smaller than 0.1 eV. 


4.3 Comparison with results of previously utilized many- 
body methods 

The foregoing analysis drew attention that both the quantum 
chemical method and the basis set utilized are important for 
correct EA- and /P-estimations. Employing small basis sets 
is particularly tempting for truly ab initio approaches, which 
are otherwise impracticable, as they require much more RAM- 
memory, disc space, and computational time than ubiquitous 
DET-flavors. 

In this subsection, we will scrutinize the reliability of the 
results obtained by post-DET approaches reported in two ear¬ 
lier studies-22i^. To assess the validity of those methods, we 
will present a comparison with results of the present methods 
also using the same small basis sets of ref. 22 i^. 

Ref.— reported results obtained within an MP2-like ap¬ 
proach for two molecules, namely (S-CgH 4 -S and S-CgH 3 E- 
S), which are similar to BDT(=HS-CgH 4 -SH). In Table[3 we 
present results for these molecules obtained within the meth¬ 
ods described above along with those extracted from ref.—. 
Eollowing ref.—, we refer to the latter results as “MP2-based”. 
Still, for clarity, we should note that, in the present terminol¬ 
ogy, the method of ref.— coincides with that labeled “ 2 nd- 
order pole” here, as it corresponds to the second-order correc¬ 
tion in the the electron-electron interaction to the electronic 
self-energy ^^i^^ . So, the results labeled “MP2-based” and 
“2nd-order pole” in Table [ 8 ] should coincide. They should 
but they do not coincide; or, more precisely, they are substan¬ 
tially different. We cannot understand these large differences 
($.1 eV) visible in Table [3 The only thing not specified in 
ref.— is the molecular geometry utilized in the calculations. 
However, as actually expected, the results presented in Table 
[ 8 ]reveal that differences in (optimized) geometries have a con¬ 
siderable smaller impact (maybe ^ 0.2 eV). Most importantly, 
as seen in Table [3 it is the combined effect of an inaccurate 
method and a too small basis set that results in very large er¬ 
rors ('^ 3 eV) for IP’s. 

By using the same (DZ) basis sets as in ref.—. Table |3 
demonstrates that the results deduced within the GW method 
are not adequate to estimate the energies of the frontier or¬ 
bitals with the accuracy required for molecular transport stud¬ 
ies. Differences ^ 0.5 eV between the GW and EOM-CCSD 
visible in Table |3 are too large, given the fact that for this 
molecule (BDT) the HOMO energy offset directly measured 
by ultraviolet photoelectron spectroscopy (UPS) amounts to 
~ 1 eV— . Drawbacks of GVT-based transport approaches were 
previously pointed out— . 
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4.4 Spatial distribution of the frontier molecular orbitals 


As highlighted above, both the method and the basis sets em¬ 
ployed are essential to properly estimate the energy of the 
frontier orbitals of an isolated molecule. Still, however im¬ 
portant, the level energy of an isolated molecule is not the 
whole issue. In a molecular junction, the active molecule is 
linked to electrodes, which yield shifts in energy via image 
charge effects. This effect, which is well established in sur¬ 
face science^, was embodied in recent studies on molecular 
electronics in a simplified form, namely by assuming point¬ 
like molecular orbitals midway between electrodes-^f^. This 
assumption is indeed a comfortable approximation, as it can 
readily be implemented in one-shot DFTh-E transport calcula- 
tions-^i^. 

The obvious critical point here, on which ref.— has recently 
drawn attention, is to what extent is it legitimate to approxi¬ 
mate a real molecular orbital as a point charge. Within classi¬ 
cal electrostatics, the interaction energy of an electron located 
at z with the image charges created in two infinite planar elec¬ 
trodes can be exactly expressed b y^^'^^ 




-2v/(l) + r 


Z — Zs 


V 


Zt—Z 


( 3 ) 


Here, Zs,t are the positions of the image planes (which are 
slightly shifted from the real electrodes, see, e.g, ref.— and 
citations therein), d = Zt —Zs, and i/(z) =d logr(z)/c/z is the 
digamma function. For a real molecular orbital, the image- 
driven energy shift should be computed by weighting eqn Q 
with the MO-spatial density p^o = which is deter¬ 

mined by its wave function 'Fmo 

<!>=/ dzpfDiz)(j)iiz), ( 4 ) 

^Zs 

Pu?iz) = JJ dxdy\'¥Moix,y,z)\^. (5) 

Close to electrodes, e. g., z^Zs, eqn Q recovers the classical 
expression 

(6) 

for a single image plane. This demonstrates that, in cases of 
molecular orbitals with significant spatial extension, the main 
contributions to the image-driven energy shift <I> come from 
regions close to electrodes. 

Since, except for ref.—, the spatial extension of the domi¬ 
nant molecular orbitals did not received consideration in previ¬ 
ous studies, we have decided to systematically investigate this 
aspect for typical molecules of interest for nanotransport. To 
this aim, inspecting spatial densities of (completely unphysi¬ 
cal LUMO) Kohn-Sham orbitals makes little sense, and (es¬ 
pecially LUMO) HF orbitals may represent a too crude ap¬ 
proximation. Therefore, like in ref.—, we have calculated the 


natural orbital expansion of the corresponding reduced density 
matrices at the EOM-CCSD level. 

This is the most reliable approach to characterize the spa¬ 
tial distribution of the extra electron or hole in molecules with 
n-type (LUMO-mediated) or p-type (HOMO-mediated) con¬ 
duction, and we are not aware of a similar study conducted 
in conjunction with molecular transport at this level of theory. 
For all the molecules considered, by inspecting the natural or¬ 
bital expansion, we found that the extra electron or hole is al¬ 
most entirely (>97%) concentrated in a single natural orbital. 

Most importantly from the present standpoint, we found not 
even a molecule whose dominant MO reside in a very narrow 
spatial region around the center. 

Rather than being strongly localized close to the center, in 
cases of n-conduction, we found that the natural orbital of the 
extra electron is more or less uniformly spread over the whole 
molecule. This is illustrated by the examples depicted in Fig.[T] 
and|2] (We employed Gabedit— to generate the figures pre- 



Fig. 1 The almost singly occupied natural orbital corresponding to 
the anion’s extra electron of the 44BPY*^ anion (“LUMO”) 
obtained via EA-EOM-CCSD/aug-cc-pVDZ calculations is 
delocalized over the whole molecule. 

sented in this paper.) Therefore, the difference between the 
value O/oc = —e^\n2/d obtained by setting z= {zs+Zt)/2i\\ 
eqn Q, which corresponds to a point-like MO located in the 
middle of two infinite metallic plates^, is substantially differ¬ 
ent from the value <t>r deduced via eqn (IHi by using the re¬ 
alistic natural orbital density. For the 44BPY molecule, the 
values thus obtained for the image-driven LUMO shifts are 
<!>;. = —2.13eV and O/oc = — l.lbeV. A difference of about 
1 eV is a big effect. 

We have also computed spatial distributions of the natural 
orbital of the extra hole (“HOMO”) in cation species relevant 
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Fig. 2 The almost singly occupied natural orbital of the extra 
electron in the 2BDCN anion (“LUMO”) is delocalized over the 
whole molecule. Result of EA-EOM-CCSD/aug-cc-pVDZ 
calculations. 


for molecules exhibiting p-type conduction. The examples 
presented in Fig. [3 and 0] illustrate two different behaviors, 
which we found to be characteristic for HOMO distributions. 



Fig. 3 The almost singly occupied natural orbital of the extra hole 
(“HOMO”) in the BDT cation is delocalized over the whole 
molecule. Result of IP-EOM-CCSD/aug-cc-pVDZ calculations. 

On one side, we found HOMO’s like that presented in Fig.[3l 
which are similar to the LUMO’s discussed above; they are 
delocalized over the whole molecule. On the other side, we 
found homo’s strongly localized on the anchoring groups. 
Typical for this behavior are alkanemono- and dithiols, as il¬ 
lustrated in Fig. 0] 

To obtain the above value 'I>r(= —2.16eV) for 44BPY, we 
applied the cutoff procedure near electrodes described in de¬ 
tail in ref.—. Results of preliminary calculations with metal 
atoms linked to the molecules with LUMO-mediated conduc¬ 
tion considered here indicate that, like the case of ref.—, the 
spatial density LUMO (single occupied natural orbital corre¬ 
sponding to the extra electron) does not substantially pene¬ 
trate into electrodes. This is important: corrections due to 



Fig. 4 The almost singly occupied natural orbital of the extra hole 
(“HOMO”) in the hexanemonothiol cation is localized in the 
vicinity of the sulfur atom. Result of IP-EOM-CCSD/aug-cc-pVDZ 
calculations. 


image charges are not dramatically affected by the cutoff pro¬ 
cedure close to electrodes. This behavior contrasts to that of 
the homo’s. Whether delocalized (like that of Fig.O or lo¬ 
calized (like that from Fig. 01i on the terminal groups of the 
isolated molecules, we found that the HOMO distributions 
substantially penetrates into electrodes. Although we only 
checked that this happens for molecules of the types discussed 
above, we believe that this is a general HOMO property that 
ensures stable molecule-electrode bonds. A cutoff procedure 
is needed to eliminate spurious divergences of the classical 
expression of the interaction energy with image charges at 
2 = Zs,t {cf- eqn (O and (|6l)], which ignores quantum mechani¬ 
cal effects and electrodes’ atomistic structure. A possible cut¬ 
off procedure consists of multiplying the classical expression 
with factors 1 - exp(-^ \z - Zsj) (°= \z - Zs,t | as z ->■ Zs,r); see, 
e.g, ref.— and citations therein. Cutoff procedures make sense 
only if the results do not sensitively depend on the value of 
cutoff parameter /r, and this cannot be the case if (HO)MO 
densities have non-negligible values in spatial regions close to 
electrodes (z Zs,t)- 

4.5 The exciton binding energy as evidence for important 
electron correlations 

As a possible way to quantify electron correlations, the solid- 
state community employs the difference between the so-called 
charge gap and the optical gap(=lowest excitation energy) 
A„; see e.g. ref. and citations therein. The charge gap, 
which is what molecular physicists normally call the HOMO- 
LUMO gap, can be expressed as A^ = IP — EA. Loosely 
speaking [because in reality the single-particle (MO) picture 
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breaks down in the exciton problem], the difference between 
the charge gap and the optical gap is that, in determining EA, 
both HOMO and LUMO are occupied. By contrast, in the 
lowest optical excitation, the HOMO becomes empty as the 
LUMO becomes occupied. A^. should be larger than A„ be¬ 
cause of the (negative) attraction energy between the oppo¬ 
sitely charged electron (LUMO) and hole (HOMO) This dif¬ 
ference is referred to as the exciton binding energy EBE (see, 
e.g., ref.^) 

EBE = Ac-Ao. (7) 

The various tables of the present paper only include results 
for EBE obtained within the most accurate method utilized 
(EOM-CCSD). The EBis-values shown there are substantial, 
amounting to up to ^ 50% of the charge gap. As expected 
for molecules with aromatic units and delocalized electrons, 
the EBE decreases with increasing molecular size; compare 
the £B£-values for BOON (Table E]) and 2BDCN (Table [Hi, 
and for BDT (TableH and 2BDT (TablejS]). The £B£-values 
{EBE ~ 3.5 —4.6eV) estimated for all the molecules analyzed 
in the present paper are substantially larger than, e.g., for n- 
conjugated organic thin films {EBE ^ 0.6 — L4eV)^. 

So, one should conclude that, for species of interest for 
molecular electronics, electron correlations are very strong. 
This aspect may be quite relevant for developing correlated 
transport approaches^; e.g., even if the charge transport is 
dominated by the LUMO, an electron traveling through the 
molecule can interact with the HOMO. 

5 Conclusion 

In this paper, we have presented benchmark quantum chemical 
calculations for the lowest electron attachment and ionization 
energies of several isolated molecular species of interest for 
molecular electronics. In assessing the importance of the ac¬ 
curacy to be achieved by estimating EA’s and IP’s, it is worth 
mentioning that a proper understanding of the charge trans¬ 
port at nanoscale does not only mean to reproduce the (order 
of) magnitude of the currents (which can be adjusted by “ma¬ 
nipulating” both £o and the broadening functions T not consid¬ 
ered here), but also the biases V characterizing a non-Ohmic 
regime, which is basically determined by Eq {eV Sq, see^). 

The main results presented above can be summarized as fol¬ 
lows; 

(i) For all molecules, the differences between the HF-MO 
energy (Koopmans theorem) and A-SCF values are large. This 
demonstrates that orbital relaxation is substantial. Electron 
correlations are also important, as revealed by important de¬ 
partures from the SCF results as well as by substantial differ¬ 
ences between the various post-SCF methods considered. 

(ii) The present results demonstrate the need both for accu¬ 
rate methods and good basis sets beyond those currently uti¬ 


lized in transport approaches. In particular, employing basis 
sets with sufficient diffuse functions is essential to correctly 
describe electron affinities. 

(iii) Kohn-Sham orbital energies can by no means be used to 
estimate IP’s and EA’s. Even with A-corrections, DFT-based 
methods do not appear to achieve the desired accuracy of es¬ 
timating the relevant MO energy offsets eo. As visible in the 
various tables, because A-DFT estimates are much weakly de¬ 
pendent on the basis sets that those of EOM-CCSD; so, A-DFT 
may convey a false impression on the importance of the basis 
sets to be utilized in calculations. 

(iv) MP2-based methods appear to be completely inade¬ 
quate for describing anions. For cations, they may yield 
substantially different results; e.g., compare the deduced via 
A-MP2 and the second-order order correction to self-energy 
(“2nd-order pole) (> 1.6 eV in Table |5ll. Examples showing 
that different methods to include second-order terms in the 
electron-electron interaction in other contexts were discussed 
earlier—. 

(v) For the presently investigated molecules exhibiting p- 
type conduction, the OVGF method represents an excellent 
compromise in terms of computational effort and accuracy of 
/P-estimates. Unlike the other diagrammatic methods con¬ 
sidered here {GW and ADC), the OVGF method does not re¬ 
quire to self-consistently solve a(n integral) Dyson equation; 
the electron self-energy E(e) can be expressed in closed ana¬ 
lytical form, and what needs to solve is a nonlinear algebraic 
equation for e— . To be fair, let us also mention that, for the 
presently considered molecules that form junctions character¬ 
ized by n-type conduction, the OVGF method turned out to be 
totally inadequate. 

(vi) The spatial distribution of the frontier orbitals plays an 
important role to reliably estimate image-driven shifts of the 
relevant MO-energies. To the best of our knowledge, this is 
the first systematic study on the spatial distribution of the extra 
electron or hole in molecular species of interest for molecular 
transport at this (EOM-CCSD natural orbital expansion) level 
of theory. None of the molecules considered in this paper was 
found to possess point-like frontier molecular orbitals, a fact 
that contradicts the common assumption made in the field. 
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44BPY/Method 


Basis set 


No. functions EA (eV) EBE (eV) 


DZ 

6-31G** 

svp 

cc-pVDZ 


136 

208 

208 

208 


-0.684 

- 0.686 

-0.421 


EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD 

EOM-CCSD(T) 

EOM-CC2 

ADC(2) 

A-SCE 

A-MP2 

A-CCSD 

A-DPT/B3LYP 

A-DET 

A-DET/B3LYP 

A-DET 

A-DET 

A-DET 

A-DET 

A-DET 

A-DET 

A-DET/B3LYP 

A-DET 

A-DET 

-KS-LUMO 
Koopmans theorem 
2nd-order pole (98.3%) 
3rd-order pole (98.6%) 
OVGF (98.5% ) 


dzp 

220 

6-31IG** 

264 

tzp 

276 

TZ2P 

360 

qz2p 

416 

cc-pVTZ 

472 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

DZ 

136 

6-3IG* 

196 

svp 

208 

6-31G** 

208 

cc-pVDZ 

208 

6-31-i-G(d) 

244 

6-311G** 

264 

6-31-i-G(d,p) 

268 

6-31-i-i-G(d,p) 

276 

cc-pVTZ 

472 

aug-cc-pVDZ 

348 

aug-cc-pVTZ 

736 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 

aug-cc-pVDZ 

348 


-0.484 

-0.353 

-0.340 

-0.290 

-0.0247 

-0.0123 

-0.438 

0.0322 4.487 

0.0293 

0.360 

0.370 

-0.564 

-0.500 

0.0043 

0.185 

-0.126 

0.124 

-0.115 

0.0290 

0.396 

0.188 

0.407 

0.407 

0.281 

0.444 

0.467 

2.068 

-0.832 

-0.466 

-0.562 

-0.540 


Table 1 The lowest electron affinity {EA) of 4,4'-bipyridine (44BPY) computed by various methods and the exciton binding energy {EBE). 
The geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 







BDCN/Method 

Basis set 

No. functions 

EA (eV) 

EBE (eV) 

EOM-CCSD 

DZ 

108 

-0.0833 


EOM-CCSD 

6-31G** 

170 

0.0893 


EOM-CCSD 

aug-cc-pVDZ 

266 

0.717 

4.614 

EOM-CCSD(T) 

aug-cc-pVDZ 

266 

0.717 


EOM-CC2 

aug-cc-pVDZ 

266 

1.047 


ADC(2) 

aug-cc-pVDZ 

266 

1.107 


A-MP2 

aug-cc-pVDZ 

266 

0.394 


A-CC2 

aug-cc-pVDZ 

266 

0.602 


A-CCSD 

aug-cc-pVDZ 

266 

0.678 


A-DET/B3LYP 

aug-cc-pVDZ 

266 

1.127 


A-DET/B3LYP 

aug-cc-pVTZ 

552 

1.155 


A-SCE 

aug-cc-pVDZ 

266 

0.350 


-KS-LUMO 

aug-cc-pVDZ 

266 

2.918 


Koopmans theorem 

aug-cc-pVDZ 

266 

-0.536 


2nd-order pole (88.6%) 

aug-cc-pVDZ 

266 

1.034 


3rd-order pole (91.1%) 

aug-cc-pVDZ 

266 

0.253 


OVGF (90.6% ) 

aug-cc-pVDZ 

266 

0.445 



Table 2 The lowest electron affinity (£A) of 1,4-dicyanobenzene (BDCN) computed by various methods and the exciton binding energy 
(EBE). The geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 


2BDCN/Method 

Basis set 

No. functions 

EA (eV) 

EBE (eV) 

EOM-CCSD 

DZ 

176 

-0.113 


EOM-CCSD 

6-31G** 

264 

0.0677 


EOM-CCSD 

aug-cc-pVDZ 

440 

0.697 

3.567 

EOM-CCSD(T) 

aug-cc-pVDZ 

440 

0.697 


EOM-CC2 

aug-cc-pVDZ 

440 

1.033 


A-MP2 

aug-cc-pVDZ 

440 

0.0227 


A-CCSD 

aug-cc-pVDZ 

440 

0.601 


A-CC2 

aug-cc-pVDZ 

440 

0.611 


A-DET/B3LYP 

6-311-H-HG(d,p) 

408 

1.228 


A-SCE 

aug-cc-pVDZ 

440 

0.121 


A-SCE 

6-311-l-tG(d,p) 

408 

0.117 


-KS-LUMO 

aug-cc-pVDZ 

440 

2.685 


Koopmans theorem 

aug-cc-pVDZ 

440 

-0.617 


2nd-order pole (97.9%) 

aug-cc-pVDZ 

440 

-0.171 


3rd-order pole (98.3%) 

aug-cc-pVDZ 

440 

-0.330 


OVGF (98.2% ) 

aug-cc-pVDZ 

440 

-0.292 



Table 3 The lowest electron affinity (EA) of 4,4’-dicyano-1,1’-biphenyl (2BDCN) computed by various methods and the exciton binding 
energy {EBE). The geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 
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BDT/Method 

Basis set 

No. functions 

/P(eV) 

EOM-CCSD 

DZ 

108 

7.492 

EOM-CCSD 

6-31G** 

150 

7.669 

EOM-CCSD 

svp 

150 

7.671 

EOM-CCSD 

cc-pVDZ 

150 

7.642 

EOM-CCSD 

dzp 

166 

7.632 

EOM-CCSD 

cc-pVTZ 

332 

7.922 

EOM-CCSD 

aug-cc-pVDZ 

246 

7.816 

EOM-CCSD(T) 

aug-cc-pVDZ 

246 

7.816 

EOM-CCSD(2) 

aug-cc-pVDZ 

246 

7.883 

EOM-CC2 

aug-cc-pVDZ 

246 

7.497 

A-MP2 

aug-cc-pVDZ 

246 

8.114 

A-CCSD 

6-31G** 

150 

7.590 

A-CCSD 

aug-c-pVDZ 

246 

7.858 

A-CC2 

aug-cc-pVDZ 

246 

7.874 

A-DFT/B3LYP 

DZ 

108 

7.690 

A-DFT/B3LYP 

6-3IG* 

140 

7.527 

A-DFT/B3LYP 

svp 

150 

7.559 

A-DFT/B3LYP 

cc-pVDZ 

150 

7.549 

A-DFT/B3LYP 

6-31G** 

158 

7.531 

A-DFT/B3LYP 

6-31IG* 

196 

7.659 

A-DFT/B3LYP 

6-31-H-HG(d,p) 

196 

7.750 

A-DFT/B3LYP 

6-311-i-G(d,p) 

228 

7.759 

A-DFT/B3LYP 

6-311-i-i-G(d,p) 

234 

7.770 

A-DFT/B3LYP 

cc-pVTZ 

332 

7.640 

A-DFT/B3LYP 

aug-cc-pVDZ 

246 

7.669 

A-DFT/B3LYP 

aug-cc-pVTZ 

514 

7.671 

A-SCF 

aug-cc-pVDZ 

246 

7.149 

-KS-HOMO 

aug-cc-pVDZ 

246 

5.865 

Koopmans theorem 

aug-cc-pVDZ 

246 

8.038 

2nd-order pole (87.4%) 

aug-cc-pVDZ 

246 

7.457 

3rd-order pole (90.1%) 

aug-cc-pVDZ 

246 

7.960 

OVGF (89.3%) 

aug-cc-pVDZ 

246 

7.750 


Table 4 The lowest ionization energy {IP) of benzenedithiol (BDT) computed by various methods and the exciton binding energy (EBE). The 
geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 
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2BDT/Method 

Basis set 

No. functions 

/P(eV) 

EOM-CCSD 

DZ 

176 

7.294 

EOM-CCSD 

6-31G** 

254 

7.340 

EOM-CCSD 

svp 

254 

7.483 

EOM-CCSD 

dzp 

276 

7.437 

EOM-CCSD 

cc-pVDZ 

254 

7.433 

EOM-CCSD 

aug-cc-pVDZ 

420 

7.593 

EOM-CCSD(T) 

aug-cc-pVDZ 

420 

7.593 

EOM-CCSD(2) 

aug-cc-pVDZ 

420 

7.747 

EOM-CC2 

aug-cc-pVDZ 

420 

7.303 

A-MP2 

6-31G* 

254 

8.629 

A-MP2 

6-311-l-l-G(d,p) 

394 

8.831 

A-CCSD 

6-3IG* 

254 

7.436 

A-CCSD 

aug-cc-pVDZ 

420 

7.714 

A-CC2 

aug-cc-pVDZ 

420 

7.888 

A-DPT/B3LYP 

DZ 

176 

7.255 

A-DPT/B3LYP 

6-3IG* 

238 

7.092 

A-DPT/B3LYP 

svp 

254 

7.172 

A-DPT/B3LYP 

cc-pVDZ 

254 

7.149 

A-DPT/B3LYP 

6-31G** 

268 

7.097 

A-DPT/B3LYP 

6-31IG* 

298 

7.248 

A-DPT/B3LYP 

6-31IG** 

328 

7.256 

A-DPT/B3LYP 

cc-pVTZ 

568 

7.242 

A-DPT/B3LYP 

aug-cc-pVDZ 

420 

7.269 

A-SCF 

aug-cc-pVDZ 

420 

6.671 

A-SCF 

6311-l-hG(d,p) 

394 

6.819 

-KS-HOMO 

6-311-i-i-G(d,p) 

394 

5.825 

Koopmans theorem 

6-311-i-i-G(d,p) 

394 

7.778 

2nd-order pole (86.5%) 

6-311-i-i-G(d,p) 

394 

7.260 

3rd-order pole (89.8%) 

6-311-i-i-G(d,p) 

394 

7.629 

OVGF (88.8% ) 

6-311-i-i-G(d,p) 

394 

7.522 


Table 5 The lowest ionization energy (IP) of dibenzenedithiol (2BDT) computed by various methods and the exciton binding energy (EBE). 
The geometry of the neutral molecule has been optimized at DFT/B3LYP/6-31++G(d,p) level—. 
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C6MT/Method 

Basis set 

No. functions 

/P(eV) 

EBE (eV) 

EOM-CCSD 

DZ 

106 

8.543 


EOM-CCSD 

6-31G** 

172 

8.741 


EOM-CCSD 

svp 

172 

8.806 


EOM-CCSD 

dzp 

183 

8.785 


EOM-CCSD 

cc-pVDZ 

172 

8.803 


EOM-CCSD 

cc-pVTZ 

291 

9.107 


EOM-CCSD 

aug-cc-pVDZ 

291 

8.966 

4.130 

EOM-CCSD(T) 

aug-cc-pVDZ 

291 

8.966 


EOM-CCSD(2) 

aug-cc-pVDZ 

291 

8.799 


EOM-CC2 

aug-cc-pVDZ 

291 

8.806 


ADC(2) 

aug-cc-pVDZ 

291 

8.623 


A-MP2 

aug-cc-pVDZ 

291 

9.0448 


A-CCSD 

aug-cc-pVDZ 

291 

8.922 


A-CC2 

aug-cc-pVDZ 

291 

9.012 


A-DPT/B3LYP 

DZ 

106 

8.989 


A-DPT/B3LYP 

svp 

172 

8.925 


A-DPT/B3LYP 

6-31G* 

137 

8.930 


A-DPT/B3LYP 

cc-pVDZ 

172 

8.918 


A-DPT/B3LYP 

6-31G** 

179 

8.926 


A-DPT/B3LYP 

6-31IG* 

176 

9.019 


A-DPT/B3LYP 

cc-pVTZ 

410 

8.975 


A-DPT/B3LYP 

aug-cc-pVDZ 

291 

8.994 


A-DPT/B3LYP 

aug-cc-pVTZ 

648 

8.990 


A-SCE 

aug-cc-pVDZ 

291 

7.988 


-KS-HOMO 

aug-cc-pVDZ 

291 

6.490 


Koopmans theorem 

aug-cc-pVDZ 

291 

9.575 


2nd-order pole (90.6%) 

aug-cc-pVDZ 

291 

8.632 


3rd-order pole (91.1%) 

aug-cc-pVDZ 

291 

9.109 


OVGF (90.7% ) 

aug-cc-pVDZ 

291 

9.058 



Table 6 The lowest ionization energy {IP) of hexanemonothiol (C6MT) computed by various methods and the exciton binding energy (EBE). 
The geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 
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C6DT/Method 

Basis set 

No. functions 

/P(eV) 

EOM-CCSD 

DZ 

124 

8.591 

EOM-CCSD 

6-31G* 

148 

8.748 

EOM-CCSD 

svp 

190 

8.847 

EOM-CCSD 

dzp 

206 

8.829 

EOM-CCSD 

cc-pVDZ 

190 

8.846 

EOM-CCSD 

cc-pVTZ 

444 

9.162 

EOM-CCSD 

aug-cc-pVDZ 

318 

9.026 

EOM-CCSD(T) 

aug-cc-pVDZ 

318 

9.026 

EOM-CCSD(2) 

aug-cc-pVDZ 

318 

8.859 

EOM-CC2 

aug-cc-pVDZ 

318 

8.482 

ADC(2) 

aug-cc-pVDZ 

318 

8.696 

A-MP2 

aug-cc-pVDZ 

318 

9.092 

A-CCSD 

aug-cc-pVDZ 

318 

8.974 

A-CC2 

aug-cc-pVDZ 

318 

9.063 

A-DPT/B3LYP 

DZ 

124 

8.209 

A-DPT/B3LYP 

svp 

190 

8.111 

A-DPT/B3LYP 

cc-pVDZ 

190 

8.136 

A-DPT/B3LYP 

6-31G** 

198 

8.167 

A-DPT/B3LYP 

6-31IG** 

244 

8.287 

A-DPT/B3LYP 

cc-pVTZ 

444 

8.246 

A-DPT/B3LYP 

aug-cc-pVDZ 

318 

8.273 

A-DPT/B3LYP 

aug-cc-pVTZ 

698 

8.279 

A-SCE 

aug-cc-pVDZ 

318 

8.042 

-KS-HOMO 

aug-cc-pVDZ 

318 

6.535 

Koopmans theorem 

aug-cc-pVDZ 

318 

9.628 

2nd-order pole (90.7%) 

aug-cc-pVDZ 

318 

8.712 

3rd-order pole (91.1%) 

aug-cc-pVDZ 

318 

9.158 

OVGF (90.7% ) 

aug-cc-pVDZ 

318 

9.116 


Table 7 The lowest ionization energy (IP) of hexanedithiol (C6DT) computed by various methods and the exciton binding energy (EBE). The 
geometry of the neutral molecule has been optimized at DFT/B3LYP/aug-cc-pVDZ level. 
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Molecule/Optimization 

Method 

Basis set 

No. functions 

IP(eV) 

BDT 

optimized B3LYP/aug-cc-pVDZ 

2nd-order pole (87.4%) 

aug-cc-pVDZ 

246 

7.457 


3rd-order pole (90.1 %) 

aug-cc-pVDZ 

246 

7.960 


OVGF (89.3%) 

aug-cc-pVDZ 

246 

7.750 


2nd-order pole (91.9%) 

STO-6G 

54 

4.702 


3rd-order pole (91.5%) 

STO-6G 

54 

4.775 


OVGF (91.4%) 

STO-6G 

54 

4.754 

S-CgH.-S 

Refi^ 

MP2-based 

STO-6G 

52 

6.5267 

optimized B3LYP/aug-cc-pVDZ 

FOM-CCSD 

aug-cc-pVDZ 

228 

8.920 


2nd-order pole (87.3%) 

STO-6G 

52 

5.139 


3rd-order pole (83.6%) 

STO-6G 

52 

5.319 


OVGF (82.9%) 

STO-6G 

52 

5.353 

optimized MP2/STO-6G 

2nd-order pole (88.3%) 

STO-6G 

52 

5.382 


3rd-order pole (85.0%) 

STO-6G 

52 

5.499 


OVGF (84.4%) 

STO-6G 

52 

5.521 

S-CgH^F-S 

Ref.^ 

MP2-based 

STO-6G 

56 

6.8019 

optimized B3LYP/aug-cc-pVDZ 

FOM-CCSD 

aug-cc-pVDZ 

242 

8.229 


2nd-order pole (87.4%) 

STO-6G 

56 

5.235 


3rd-order pole (83.7%) 

STO-6G 

56 

5.370 


OVGF (83.0%) 

STO-6G 

56 

5.395 

optimized MP2/STO-6G 

2nd-order pole (88.2%) 

STO-6G 

56 

5.444 


3rd-order pole (85.1%) 

STO-6G 

56 

5.522 


OVGF (84.5%) 

STO-6G 

56 

5.534 


Table 8 MP2-based results for the lowest ionization energy (IP) of the benzenedithiol (BDT) and BDT-like molecules considered in ref.— 
along with results obtained within methods utilized in this paper. 


Method 

IP(eV) 

EA (eV) 

HOMO-LUMO gap (eV) 

FOM-CCSD 

7.492 

-2.016 

9.500 

OVGF 

7.560 

-2.073 

9.633 

GW 

6.9 

-2.2 

9.1 


Table 9 Results for the lowest ionization energy (IP), electron affinity (EA), and HOMO-LUMO gap obtained by using several methods and 
the same double zeta (DZ) basis sets. Notice that the GW estimates taken from ref.— significantly differ from the EOM-CCSD values, which 
are very accurately estimated within the OVGF. One should remark that the fact that all methods predict a stable anion (EA > 0) is an artefact 
of the small (DZ) basis set utilized in this table. 
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